Epicyclic oscillations of non-slender fluid tori 
around Kerr black holes 



Odele Straub 1 and Eva Sramkova 2 

1 Copernicus Astronomical Centre PAN, Bartycka 18, 00-716 Warsaw, Poland 
2 Department of Physics, Silesian University in Opava, Bezrucovo nam. 13, 746-01 
Opava, Czech Republic 

E-mail: odele@camk.edu.pl, sram_eva@centrum.cz 

Abstract. Considering epicyclic oscillations of pressure-supported perfect fluid tori 
orbiting Kerr black holes we examine non-geodesic (pressure) effects on the epicyclic 
modes properties. Using a perturbation method we derive fully general relativistic 
formulas for eigenfunctions and eigenfrequencies of the radial and vertical epicyclic 
modes of a slightly non-slender, constant specific angular momentum torus up to 
second-order accuracy with respect to the torus thickness. The behaviour of the 
axisymmetric and lowest-order (m = ±1) non-axisymmetric epicyclic modes is 
investigated. For an arbitrary black hole spin we find that, in comparison with 
the (axisymmetric) epicyclic frequencies of free test particles, non-slender tori receive 
negative pressure corrections and exhibit thus lower frequencies. Our findings are in 
qualitative agreement with the results of a recent pseudo- Newtonian study of analogous 
problem defined within the Paczyhski-Wiita potential. Implications of our results on 
the high-frequency QPO models dealing with epicyclic oscillations are addressed. 
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1. Introduction 



Oscillations of black hole accretion discs have been studied extensively in various 
astrophysical contexts. Investigations of hydrodynamic oscillation modes of 
geometrically thin accretion discs revealed three fundamental, discoseismic classes 
of modes: acous ti c pressure p-mo d es, gr a vity g- modes and corrugati o n c-m odes. 
Kato et al I Jl998h : IWaeonerl Jl999h : iKatol fhoOlah and IWaeoner et al\ fhoOlh give 
comprehensive reviews on the subject of relativistic 'discoseismology'. Geometrically 
thick discs (tori) have been examined rather less and mainly in matters of their stability 



(e.g. 



Papaloizou fc Pringld . I1984J ; iKojimal . Il986l ). Only recently oscill atory modes of 



fluid tori were explored in more detail in several numerical studies (e.g. IZanotti et al 
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Montero etal 



200 31; IRezzolla et all l2003al : iRubio-Herrera and Led . l2005alJbl; ISramkova et al 



2007; 



20071 ) and in a purely analytic work by lBlaes. Arras and Fragile! ( 120061 ). 



The latter presents a thorough analysis of oscillatory modes of relativistic slender tori. 

The analysis of disc oscillation modes is motivated by observations of quasi- 
peri odic oscillations (QPOs) in the l i ght curves of Galac tic low-mass X-ray binaries 



McClintock and Remillard. 



2003; 



van der Klisl . I2Q04J . for reviews). In particular 



(see 

the understanding of high frequency (HF) QPOs, which are presumably a strong gravity 
phenomenon, would provide a deeper insight into the innermost regions of accretion discs 
and the very nature of compact objects. Black hole HF QPOs occur at frequencies that 
are constant in time and characteristic for a particular source. If more than one HF QPO 
is detected in a given system, the frequencies typically appe ar in ratios of small natural 
numbers, whereas in most cases the r atio is close to 3:2 (lAbramowicz and Kluzniakl . 



200ll : iMcClintock and Remillardl . 120031 ). 

Several models ex plain HF QPOs in terms of disc (or blob) oscillations. 
Stella and Vietril (119981 ) for instance consider an orbiting hot spot and propose that 



HF QPOs arise due to a modulation of the spot radiation by its precessional motion. 
Pointi ng out the observed rational ratios of HF QPO pairs iKluzniak and Abramowicz 
( 120001 ) and lAbramowicz and Kluzniakl (120011 ) suggest an underlying non-linear 
resonance between some modes of accretion disc oscillation. In the resonance model the 
mode pair is commonly represented by the radial and vertical epicy clic oscillation. The 
importance of epicyclic oscillations is also stressed by iKatol (l2001bl ) who attributes the 
origin of HF QPOs to non-axisymmetric g-mode oscillations. The corotation resonance 
in a disc which is deformed by a warp would be responsible for the excitation of g- 
modes and general relativistic effects trap them near the inner edge of the accretion disc 
( IKatol . 



2003 



2004 ) . HF QPOs could also result from acoustic p-mode oscill ations of a 



small accretion torus orb i ting c lose to the black hole (IRezzolla et al I . l2003bl ). Recently 



Blaes. Arras and Fragile! (120061 ) discussed the possibility that the vertical epicyclic and 



the lowest-order (acoustic) breathing mode of a relativistic slender torus might represent 
the two black hole HF QPOs in 3:2 ratio. 

The main interest of our work lies in the epicyclic modes of oscillations. Most 
models that are dealing with them consider geodesic flows and are based on free test 
particles. However, non-geodesic effects related to e.g. magnetic fields, viscosity or 
pressure forces may play a certain role in this concern. The aim of this paper is to 
investigate such non-geodesic effects on the two epicyclic modes by means of a pressure- 
supported perfect fluid torus and to find the consequential pressure corrections to the 
mode eigenfunctions and eigenfrequencies. 

For an infinitely slender torus, the frequencies of epicyclic oscillations are consistent 
with the epicyclic frequencies of free test particles on a circular or bit in the equatorial 



plane , where the torus pressure exhibits a maximum value ( lAbramowicz et al 



20061 ). It h as been shown th at the epicyclic modes ma y be retained a lso for 



Rubio-Herrera and Led . 



thicker tori (jBlaes et al I. 12007). Numerical simul ations (IRezzolla et al I . 12003a 



2005a 



bj; ISramkova et al 1 . 120071 ) as well as analytic calculations 
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in pseudo-Newtonian approximation (IBlaes et al 1 . 120071 ) showed that with growing torus 
thickness the (axisymmetric) epicyclic frequencies decrease. What we have in mind here 
is a slightly non-slender torus (with a radial extent that is very small in comparison to 
its distance from the central object) in hydrostatic equilibrium in Kerr spacetime. We 
perform a second-order perturbation analysis of the eigenf unctions and eigenfrequencies 
of both modes with respect to the torus thickness and derive exact analytic formulas 
for t he pressure c orrect ions. 



Blaes et al I (120071 ) have studied an analogous pro blem in the framework o f 



Newtonian physics using the pseudo- Newtonian potential of lPaczyhski and Wiital (119801 ) 
to model the relevant general relativistic ef fects. Our w o rk rep resents a generalisation 
of their results into Kerr geometry. Like IBlaes et al I (120071 ) we assume a non-self- 
gravitating, non-magnetic, stationary torus with a constant specific angular momentum 
distribution. We neglect self-gravity because the mass of a black hole exceeds the mass 
of a torus many times over. Although effects of magnetic fields are neither negligible 
nor irrelevant to our problem, we study here a purely hydrodynamical case. Our 
calculations are only valid as long as magnetic fields are unimportant, i.e., as long as the 
torus pressure dominates the magnetic pressure. The stability of these hydrodynamic 
modes in presence of magnetohydrodynamical turbulence is an issue that needs yet 
to be investigated. The above assumptions are supported by numerical simulations. 



Proga fc Begelmanl (120031 ) show for instance in an inviscid hydrodynamical simulation 



that the inner axisymmetric accretion flow settles into a pressure-rotation supported 
torus with constant specific angular momentum. Once magnetic fields are introduced 
the angular momentum distribution gets a different profile, torus-l ike configurations are , 



however, still seen a s "inner tori" in global MHD simulations (e.g. iMachida et al 1 . 12006 



Fragile et al 1 . 120081 ) . Although the described torus setup is not the most likely to be 



found in nature, we think it is reasonable to assume such a configuration as a first 
approximation. 

In this work we focus on a mathematical description of the problem, whereas the 
astrophysical applications shall be presented separately. The paper is outlined in the 
following manner: In sections [2] and [3] we give a brief introduction to the problem writing 
down the equations that describe the relativistic equilibrium tori and the relativistic 
Papaloizou-Pringle equation. Then, in section HI we describe the perturbation method 
and derive formulas for the radial and vertical epicyclic mode eigenfunctions and 
eigenfrequencies. The results are presented for different values of the black hole spin 
parameter in section [5] and discussed in section [6j 



2. Equilibrium configuration 

Consider an axisymmetric, non-self-gravitating perfect fluid torus in hydrostatic 
equilibrium on the background of the Kerr geometry. The flow of fluid is stationary 
and in a state of pure rotation. Generally, the line element of a stationary, axially 
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symmetric spacetime is given in Boyer-Lindquist coordinates (t, r, 6, 0) by 

ds 2 = g tt dt 2 + 2g ttj) dtd<f> + g„dr 2 + g 9e d9 2 + g^dcj) 2 . (1) 

We take the ( h ++) signature and units where c = G = M = 1. The explicit 

expressions for the covariant and contravariant coefficients of the Kerr metric then write 

9u = ~ (l - 9 tt 



9t<f> 


E 




E 


9rr 


~ A' 




= s, 






9<j>4> 


= (V 2 



Sy EA' 

sin^, q t<p = -, 

y EA 

r = f , 

9 2a 2 r . 9 A . o „ A - a 2 sin 2 9 

a + — »)-»•«. sw fl ■ P) 

where E = r 2 + a 2 cos 2 0, A = r 2 — 2r + a 2 , H = (r 2 + a 2 ) 2 — a 2 A sin 2 6 and M is the 
mass and a the specific angular momentum (spin) of the black hole. 

Because the flow is assumed to be purely azimuthal, the four-velocity has only 
two non-zero components, m m = (it*, 0, 0, v/). One may derive the fluid specific angular 
momentum I, angular velocity Q, specific energy £ and the contravariant t-component 
of the four-velocity, often denoted as A, in the form 

, _ _ }H _ 9t<j> + ^9<j»j) /g\ 

u t g tt + ttg t(/> ' 

u<i> _ g** - Ig** 

u t g tt _ i g bt, ' ^ > 

S = -u t =(-g* + 2lg t *-l 2 g^)- 1 / 2 , (5) 

A= u l =(-g a -2n gt<f> -n 2 g H )- 1 / 2 . (6) 

The perfect fluid is characterised by the stress-energy tensor T^ u = (p + e)u^u v + 
pg^ v . We restrict our consideration to polytropic flows such that, measured in the fluid's 
rest frame, pressure p, internal energy density e and rest mass density p are related by 
p = Xp("+ 1 )/ n and e = np + p, where n is the polytropic index and K the polytropic 
constant. 

We assume the specific angular momentum to be constant throughout the torus, 
i.e., l(r, 9) = l — const. Such a configuration is governed by the relativistic Euler 
equation which may be written as 

= M (In £), fie{r,9}. (7) 

p + e 

Introducing the enthalpy H = J dp/(p + e), the integration of ((Tj) leads for a barotropic 
fluid, for which p = p(e), to the following form of the Bernoulli equation 

H + In 8 = const. (8) 
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that determines the structure of the torus in the r — 6 plane. The subscript zero refers 
to the special location r = r in the equatorial plane where the pressure gradients vanish 
(p has a maximal value) and the fluid moves along a geodesic line. 

For a small torus cross-section, when the adiabatic sound speed defined at 



one may write 



the pressure max imum c 2 
H ~ (n + 1)p//Q (lAbramowicz et al 1 . 120061 ). Following lAbramowicz et al I (120061 ) and 



(n + l)pn/ (npn) satisfi es c 2 << c 2 



7M). 



(9) 



Blaes et al I (120071 ) we introduce the function f(r, 9), which takes constant values at the 
isobaric and isodensity surfaces, by 
P _ Po 
P Po' 

Form the Bernoulli equation (jSJ) with the constant evaluated at the pressure maximum 
r and the above form of H we get 

1 



nc 



(ln£-ln£ ) 



(10) 



s0 



2.1. Epicyclic oscillations 

It is advantageous to introduce the effective potential 
U = g u - 2 + 

that has its minimum at the torus pressure maximum ro. A small perturbation of a 
test particle orbiting on a geodesic line r = r with I = l results in radial and vertical 
epicyclic oscillations around the equilibrium point r n at a radial uj t q and vertical ojq$ 
epicyclic frequency given by (e.g.. lAbramowicz et al I (120061 )) 

, 1 / S 2 d 2 U\ , 9 1 / S 2 d 2 U* 

2 \A 2 g rr dr 2 



vO 



and 



6*0 



2 \A 2 g ee d9 2 



(12) 



In Kerr geomet ry ([2 



Nowak and Lehr. 1998 



the above definitions lead to (Aliev and Galtsov. 1981 



Torok and Stuchlfki . 120051 ) 



U). 



rO 



n 2 



+ 



8a 3a 2 



r 



3/2 



and 



00 



Q 2 



1 - 



4a 3a 2 
+ 



,3/2 



(13) 



where f2 is the angular velocity at the pressure maximum r that in Kerr geometry 
reads Qq = 1/ (r^ 2 + a). 

In order to investigate t he behaviour of the equipo tential function / in close vicinity 
of the equilibrium point r , lAbramowicz et al I (120061 ) introduced local coordinates, 

r — ro\ , n ( 7r/2 — 6" 



x 



V9 



rrO 



and 



(14) 



y = at r . For small x and y, and a constant specific angular momentum 



satisfying x 

torus, the equipotential function / can be expressed as 
/ = 1 - ^2 (^o^ 2 + ^eoV 2 ) , 



(15) 
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where u r0 = u r0 /Q , u eQ = u eQ /Q , and (3 is a dimensionless parameter given by 

?nr 2 

P 2 - (16) 



which determines the thickness of the torus. 

Equation (Tl5l) describes the equipotential function in the vicinity of r in terms 
of the test particle epicyclic frequencie s and is congruent with the formula derived in 
Newtonian theory (see equation (9) in iBlaes et al I (120071 )). In the slender torus limit 
(3 — > 0, the torus reduces to an infmitesimally slender ring at the pressure maximum r . 
In a Newtonian 1/r potential, u r o = u eo = Q and the slender torus cross-section has 
a circular shape. In the general case, however, cu r0 ^ ojqq and the isobaric surfaces are 
ellipses with semi-axes being in the ratio of the two epicyclic frequencies. 



3. Perturbation equation 



We consider small linear perturbations around the axisymmetric and stationary torus 
equilibrium with azimuthal and time dependence in the form oc exp[i(m0 — out)]. 
The differential equation governing such perturba tions for constant specific a ngular 
momentum tori was in Newtonian theory derived by 
it was expressed in terms of a scalar variable W. 

derived the general relativistic form of the Papaloizou-Pringle equation in terms of 

dp 



Papaloizou &: Pringld (119841) where 



Abramowicz et al 



W = — 

Ap(uj — mVt) ' 

which is related to the Eulerian perturbation in the four-velocity as 



(120061 ) recently 
(17) 



Su, 



-d^W, /i€{r,0}. 



p + e 

The relativistic Papaloizou-Pringle equation writes 
1 



;is) 



-g)l/2 



[(-g) 1/2 g^f n d v W] } - {m 2 g^ - 2mug^ + u 2 g tt ) f n W 



2nA{u - mVtf 



t~ v W (19) 



where {/i, u} G {r, 6}, A = A 2 /Aq, Q = fl/flo, u = u/Q and g denotes the determinant 
of the metric. 

Following Blaes et al I (120071 ) we write (HH as 

LW = -2nA (io-mn) 2 W } 

where L is a linear operator given by 



(20) 



L = [ d,(^g-)g^fd u + d^nfdu + g^nd,fd u + g^f&, 



(m 2 g^ - W + u 2 g tt ) f}f3 2 r 2 . 



(21) 
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4. Expanding the relativistic Papaloizou-Pringle equation about the 
slender torus limit 

We now use a perturbation method to derive the expressions for eigenfunctions and 
eigenfrequencies of the epicyclic modes for thicker tori {(3 > 0). 

We start by transforming all variables to local coordinates, x = x/(3 and y = y/(3, 
measured from the equilibrium point. Then we expand u, W, A, Q, f and L into a 
power series in (3 by writing 

Q = Q(°) + QW/3 + QW[3 2 + . . . where Q e {lj, W, A, Cl, f, L}. (22) 

Substituting all variables (1221) into the perturbation equation (120]) and comparing terms 
of the same order in (3 we obtain formulas for the respective corrections to eigenfunctions 
and eigenfrequencies of the desired modes. 



4-1. Slender torus limit 



In the slender torus limit (/3 
reduces to 



0) the relativistic Papaloizou-Pringle equation ffl9l) 



/ 



(o) 



Q2 W (0) Q2 W (0) 
+ 



dx 2 



dy 2 



n 



/ g/(Q) dW<® df<® dW<® \ 



\ dx dx 



dy J 



It represents the zeroth-order of ( 120|) and may be written in operator form, 
£(0)^(0) = -2nA i0) a 2 W {0 \ 



(23) 



(24) 



where a = to^ ' — mQ^ ' denotes the zeroth-order eigenfrequency in the corotating frame 
scaled with the orbital velocity Qq. The encountered zeroth-order expansion terms are 

^ (0) = 1, (25) 
fi (0) = 1, (26) 

/ (0) = 1 - ufa? - ^eof, (27) 



d 



dx z 



d 2 df<® d dfW d 



+ ri- 



dy 2 dx dx ' dy dy ^ 
Equations (|23|) and (I24"I) are of the same form as the Newtonian slender tor us limit 
of the Papaloizou-Pringle equation (lAbramowicz et al 1 . 120061 ; iBlaes et al 1 . 120071 ). They 
represent an eigenvalue equation for a with being a self-adjoint operator with respect 
to the inner product 




(fi^W^wPdxdy^S^, 



(29) 



where the integrals are taken over the slender torus cross-section where > 0. This 
implies that the eigenvalues a are real numbers and the zeroth-order eigenfunctions 
W(o) form a complete orthonormal set. Therefore, any function defined over the torus 
cross-section may be expanded in terms of W^°\ 
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The explicit expressions for eigenfunctions and eigenfrequencies of the complete 
set of normal modes of a constant specific angular mo mentum slend er torus in 



Newtonian potential, where uj r n = ujqq = 1, were given in iBlaed (119851 ). Recently, 



Blaes. Arras and Fragile! (120061 ) derived the eigenfunctions and eigenfrequencies of the 
lowest-order modes of a general relativistic slender torus, where a> r0 ^ ujqo, for arbitrary 
specific angular momentum distribution. The eigenfunctions and eigenfrequencies of the 
simplest modes of the relativistic torus with constant specific angular momentum which 
are relevant to our calculations are specified in table [U 

In the slender torus limit the two epicyclic modes correspond to % = 1 and i = 2 
from table [U The corresponding eigenfunctions take the form 

= Wl 0) = axxe^-^ and = = a 2 ye* n +- u ™ t \ (30) 

They describe the global, incompressible modes in which the entire torus moves in 
purely radial (W^) or purely vertical (W^) direction at frequencies which are, in 
the corotating frame, consistent with the epicyclic frequencies of free test particles, 
^Vo = o"i^o = — rntt^ and ojqq = 02^0 — ^2°^ — mQ(°\ 

4-2. Non-slender torus 

4-2.1. First-order corrections Expanding the Papaloizou-Pringle equation ( 1201) to first- 
order in (3 we arrive at 

-2na*wU. (31) 

The perturbed basis of eigenfunctions may now be expressed in terms of the 

orthonormal zeroth-order basis as 

^ W = EM*f- (32) 
3 

The coefficients may be determined by taking the inner product of (l3~il with a zeroth- 
order eigenfunction W% . If the subscript k refers to a different mode from the one we 
are interested in (i.e., k 7^ i) we find 

(W} 0) \LV> - Anma&W + 2na 2 i A^\W^ ) ) 

kj = 2n{a) - a?) ' (33) 

Whereas if k refers to either of the two epicyclic modes (k = i) we obtain the formula 
for the first-order correction to the radial (i = 1) or vertical (i = 2) eigenfrequency, 

/ w (0) |LW - 4nma0V + 2nafA^\W l {0) ) 
= -± - A '-■ (34) 

\ We use the index i to label the modes of interest. It takes values i = 1 for the radial or i = 2 for the 
vertical epicyclic mode. 
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Table 1. Eigcnfunctions of the lowest-or der modes of a general slende r torus 
with constant specific angular momentum (jBlaes. Arras and Fragilel 2006 ). The 



normalisation constants ao — a§ are given in table [3] and the coefficients wn , W42 , 
W52 are specified in the appendix. 



i 









ao 




1 


a\x 




2 


0.2V 




3 


a 3 xy 




4 


0-4 (l - 


- wax 2 + wwy 2 ) 


5 


a 5 (l - 


- w 5 ix 2 + w 52 y 2 ) 



The first-order expansion terms of A, Q, f and L have a (x, y) dependence in the 

form 

A {1) = A n x, (35) 

= n n x, (36) 

/W = f u x 3 + faxf, (37) 
= (L wl x + L 102 x 3 + Li 03 sy 2 ) <9 2 
+ (L104X + L105X 3 + L W6 xy 2 ) <9| 

+ (L 107 + Li 08 x 2 + L W9 y 2 ) d s + L lw xydy, (38) 

where the explicit forms of the coefficients An, flu, fn, /12 an d L±ox —Lxio are specified 
in the appendix. 

Substituting ( I35l )-f l38|) int o (1331) and considering the eigenmodes derived in 



Blaes. Arras and Fragilel (120061 ) we find for the radial mode (i — 1) three non-zero 



coefficients that correspond to the modes i = 0, 4 and 5 given in table [TJ The resulting 
eigenfunction W r of the radial epicyclic mode for a slightly non-slender torus in first- 
order accuracy reads 

W r = a lX + P (C + C x x 2 + C 2 y 2 ) + 0(f3 2 ), (39) 

where C = a b 10 + a A bu + a 5 bi 5 , C\ = a^buWn + a 5 b 15 w 5 i and C 2 = a A b u w A2 + abbi^w^ 2 . 
The normalisation constants ao, a A , a 5 are listed in tableland the eigenfunction-related 
terms W41, w A2) W51, w$ 2 are specified in the appendix. The coefficients 6i , 614, &15, given 
by (|33|) . take the following forms 

610 = " 4n 2 (n + l)K-a?HV3 o ( L ^o + fr* + 2n(m - ^ r0 ) x 

x {An(m - LU r0 ) - 2mOu} + 2L 107 (™ + 1)^ 2 ] ^00} > ( 40 ) 

b lq = ^7 — — — 2 2 , 5 _ 5 [{^108 + 2nA u (m - u; r0 ) 2 

8n 2 (n + 1) (n + 2) (a 2 - <t 2 K u4 

-4m 2 nfin + 4mnO u a)ro} X {[2(n + 2)<I> 2 + 3W q i]uj 2 m 

-,2 W l,-.2 , T ftnf„ 1 oV-,2 1 TJ/ L-.2 1 9r,2 IJ/ 1 ,-,2 



+oo 2 r0 W q2 }u 2 + L 109 {[2(n + 2)lo 2 + W ql p 2 + 3oo 2 r0 W q2 }^ 



r0 
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1 


n 2 

U! r0 












2 


n 2 












3 


Q r0 + ui eo 












4 


{{2n + l){u> 2 r0 + Q 2 ) 


— [4n(n - 


-i)K 2 


-^ 2 o) 2 - 


h (u) 2 - 


h^ ) 2 ]V2} /(2n) 


5 


{(2n+l)(^ 2 + ^ 2 ) 


+ [An(n - 


-i)K 2 


-- 2 o) 2 - 




h^ ) 2 ] 1/2 }/(2n) 



Table 3. Normalisation constants of the eigenmodes in table [T] (jBlaes et al . 2007 ). 



They are calculated such that the eigenfunctions are normalised in the inner product 
(ESI). 



(ncj rQ u;g /7r) 1/2 

1 a [2(n+l)Q 2 }^ 2 

2 ao^n+l)^] 1 / 2 

3 a [4(n+l)(n + 2)w 2 w 2 ] 1 /2 

4 a {(n + 2)[a 2 - (^ 2 + Q 2 r0 )}/ [2na 2 - (2n + 1)K 2 + C 2 )] } 1/2 

5 a {(n + 2)[ ( t 5 2 - (co 2 + u, 2 r0 )}/ [2na 2 - (2n + l)(u, 2 + u, 2 )] } 1/2 



+2L 107 (n + 2)iu 2 iu 2 {[2(n + l)uo 2 r0 + W ql p 2 0o + Q 2 r0 W q2 }] , (41) 

where q = 4 or q = 5 in case of b u or 6 15 , respectively, and <t , 04, ct 4 , <7 5 are specified 
in table [21 

Similarly, for the vertical epicyclic mode (i = 2) we find one non-zero coefficient 
corresponding to j = 3 and receive the final expression for the eigenfunction of the 
vertical epicyclic mode, 

We = a 2 y + (3C 3 xy + 0((3 2 ) (42) 

with C3 = a 3 b 2 3- Again, 03 is the normalisation constant (see table [3]) and the 623 
coefficient is given by 

a 2 a 3 ir 



^23 



8n 2 (n + 2)(n + l)(af - e^)u>a,w; 



{Lno - 2n(m - cu eo )x 



x [^11(^0 - m) + 2mftn] } (43) 

with cr 2 , cr 3 listed in table [2j 

Then, using fl35l) - fl38l) . we find the terms 

AM\wl 0) ) = ai A n x 2 , „4 (1) |W 2 (0) > = a 2 A n xy, 

n^\w} 0) ) = ai Q u x 2 , Q W \W 2 (0) ) = a 2 Q n xy, 

L«|Wf )) = ai (L 107 + L W8 x 2 + L W9 f), L«|W 2 (0) ) = a 2 L lw xy. (44) 

Substituting (144)) into the formula (1341) for u\ , we obtain in the inner product for 
both modes i = 1 and i — 2 odd functions of x and such that the integration 
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over the elliptical torus cross-section yields zero. Therefore, to find the relevant 
pressure corrections to the radial and vertical mode frequencies, one needs to extend 
the expansion to second-order in torus thickness. 

4-2.2. Second-order corrections The perturbation equation ( 1201 expanded to second- 
order in (3 reads 

L«» W W + + #*>wf > = - In {a 2 W^ + (a 2 A^ 

-2m ( r i OW)W i (1) + [a 2 A {2) 

+m 2 (Q^) 2 -2ma0V +2aM 2) 

-2maiA {1) Cl {1) ]W^Y (45) 

Analogous to the first-order case one may write 

^ (2) = E C ^ (0) ( 46 ) 

3 

and take the inner product of (I4"5l) with a zeroth-order eigenfunction wjf\ For a 
subscript k referring to the mode of interest (k = i) we obtain the formula for the 
second-order correction, 

Anai 

+2na 2 A {2) ~ inmaiA^Cl^lwP' 



-(2) 



W^L® + 2nm 2 (ft«) 2 - Anma0 2 ^ 



(WPILM + 2na 2 A^ - Anmcx0^\W^ 



(0) 



(47) 



The second-order terms of A, Q, f and the L operator take now the form 

A {2) = A 21 x 2 + A 22 y 2 , (48) 

n {2) = n 21 x 2 + n 22 f, (49) 

/ (2) = /21X 4 + f 22 x 2 y 2 + f 23 y\ (50) 
L {2) = {L 201 x 4 + L 202 x 2 + L 203 x 2 y 2 + L 204 y 2 + L 205 y 4 } d\ 

+ {L 20 6^ 4 + L 2{n x 2 + L 208 x 2 y 2 + L 2m y 2 + L 2m y 4 } 9| 

+ {L 2 io^ 3 + L^nx + L 212 xy 2 } d s 

+ {L 213 x 2 y + L 214 y + L 215 y 3 } dy 

+ {L 216 + L 217 x 2 + L 218 y 2 } , (51) 

with the coefficients A21, A 22) Q 2 x, ^22, /21, ^22, ^23 and L 201 - L 218 again specified in 
the appendix. 

Inserting all first and second-order expansion terms into (1471) we gain a fairly 
long expression for Cu\ . Along with ff22l) it leads to the resulting formula for the 
eigenfrequency u r of the radial mode. In order to keep the overview we write the 
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expression in the following form 

B 2 4 

iu r = u r0 + m - j— J2 p i + °(P 3 ) ( 52 ) 



where 



Ano\ 

Pi = (W^ILP) + 2nm 2 (0«) 2 - Anma^ + 2naM (2) 

{(L212 + -^218) &r0 



4n(n + l)(n + 2)u^ 
+ [3 (L210 + £217 + 2m 2 nOy 
+2 (L211 + L 216 ) (n + 2) Q*)] co 2 9 
-Amn [3(Ai^n + ^21)^0 + ^22^0] ^Vo 

+2n (3A 21 u 2 + «4 22 ^o) ^ro }, (53) 



P 2 = b 10 (W± 0) \LW + 2naM {1) - ^m^O^ 1 ) |^ (0) ) 
(n + l)^o^ r 3 o ' 

P g _x = 6 l9 (^r|L« + 2^M (1) -4nm ( r 1 2 n«|^ )) 

= 6l9 2n(n + l)(n/ 2 T^ K -^ {" 2(n + 1) ^~ X 
x [3^o(-^io2 — -^105 + -^ios) + 2(n + 2)c^ u; 2 (L 101 — L 104 

+-^107) + ^ro(-^103 — -^106 + £l09 — £110)] + no ~g \prG [-^103 

+-^109 + 2(2 + n)cDg (Lioi + £107)] ~~ ^ro^eo [-^106 + -^110 
+2(n + 2)^ Lio4] + 3u 2 uj 2 (L W2 + £ios) — 3u)q Liq^ + 
n(a)e - ^ 2 )^ r o(^.ii^ r o - 2mQ u ) [2(n + 2)^ 2 a)g 
+3u 2 w ql + u) 2 w q2 ] } ■ (55) 

Here q = A or q = 5 in case of P 3 or P 4 , respectively. 

For the vertical epicyclic mode we find a much shorter expression, 

Q 9 = *«> + m- P 2 128n4(n - 1)2(n ((j 3 _ ^ ^ x 

x 0377 [Lno + 2nw 0o (Ai^eo - 2mO n )] {3Li 09 cD 2 + ^ 2 [L 108 
+Lno + 2L l07 {n + 2)u 2 + 2nu} eo (AiiLUg - 2mQ, 11 )] } 
- 8n 2 {n + l)(n + 2)^ cu r 3 (a 2 - a 3 2 ) {3 (L 215 + £218 ) ^ + <4) [^213 
+L217 + 2m 2 ntt 2 n + 2 (L 214 + L 216 ) (n + 2)lj 2 . ] - Amn [(.An^ii 

+Q 2 l) Ueo + 3fi 2 2^ro] ^80 + 2n (A 21 U) 2 Q + 3-4 2 2^o) ^90 } 

+ 0((3 3 ). (56) 
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5. The properties of epicyclic modes for non-slender tori 

The equations (1391) . (H2l) . (1521) and (1561) represent the epicyclic mode eigenfunctions and 
eigenfrequencies of non-slender tori as functions of r , a, m, n and (3. In this section 
we use these formulas to illustrate the behaviour for the axisymmetric (m = 0) and 
lowest-order non- axisymmetric (m = ±1) epicyclic modes for variable torus thickness 
and varying black hole spin value. We take the polytropic index n = 3 that refers to a 
radiation-pressure dominated torus % 

All figures that are related to frequency behaviour (jTK3l [MB display the frequencies 
defined as v = oj/2tt [jj. The horizontal axis in all these figures starts at the radius of 
the marginally stable orbit calculated for the appropriate black hole spin. 

5.1. Axisymmetric epicyclic modes 

In the axisymmetric (m = 0) case the ^-components vanish from the equations and 
the problem becomes symmetric with respect to the rotational axis (6 = 0). The 
corresponding axisymmetric radial and vertical epicyclic frequencies are shown in figures 
[U [2]and[3j The (3 = line always refers to the epicyclic frequency of a test particle, while 
the (5 > lines illustrate the behaviour of the two frequencies when the torus becomes 
thicker. As a result of the topology of the equipotential surfaces there is an upper limit 
on f3 for a given torus pressure maximum above which no closed equipotential surfaces 
and consequently no equilibrium tori may exist. This limit is incorporated into the 
figures as a dash-dotted line that defines the region of 'allowed frequencies' above it and 
to its right (inside the shaded region). 

For any black hole spin both axisymmetric frequencies decrease with increasing 
torus thickness. A comparison of the left (a = 0) and right (a = 0.8) panels of figures 
[H and [2] and the respective panel of figure [3] (a = 0.999) illustrates the influence of the 
black hole rotation. 

The radial frequency qualitatively retains for all spin values the same profile albeit 
the torus thickness. The radius where it becomes zero moves, however, away from the 
central object. It is interesting to note that for all values of a the frequency maximum 
for a fixed (3 is reduced by almost exactly the same relative amountqj. Moreover, for 
a non-extremely rotating black hole with a < 0.96, this is more or less true for any 
location of the pressure maximum r , while for a > 0.96 the frequencies at small radii 
tend to crowd together, as may be seen in figure [3l 

The vertical frequency for a Schwarzschild black hole changes its profile with 
increasing torus thickness from a monotonic function in the case of test particle frequency 

§ Note that varying n makes no relevant difference to the results. 

|| As it is commonly used, throughout the paper we refer to quantities u> as to 'frequencies' as well, 
although in exact notation they should be called 'angular velocities'. 

% This does not apply to the actual thickness of the torus, since a given fi implies for different a a 
different extent of the torus. 
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Figure 1. The axisymmetric radial epicyclic frequency as a function of the torus 
pressure maximum r for tori of various thickness that orbit a black hole of mass 
M = 1OM . The allowed frequencies corresponding to equilibrium tori lie inside the 
shaded region. Left In the Schwarzschild case a = 0. Right For a rotating black hole 
of a = 0.8. 
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Figure 2. The same as figure [TJ but for the axisymmetric vertical epicyclic mode. 



to a function exhibiting a maximum value (see the left panel of figureE]). For a Kerr black 
hole the non-monotonicity is already present for the test particle frequency (although 
for low to moderate spin values the frequency maximum is located at radii inside the 
marginally stable orbit) and it remains also for a thicker torus (see the right panels of 
figures [2] and [3]) . At very high spin values (a > 0.96), the frequency shape is modified in 
such a way that the frequency maxima for all tori almost coincide with the f3 = curve 
(figure [3], right panel). 

In order to illustrate the characteristic features of the flow for the modes of a thicker 
torus we use equations f|T8l) . (139]) and fT42|) to plot the corresponding poloidal velocity 
fields. The axisymmetric radial mode shows a relatively coherent flow with only slight 
deviations from the radial motion in the outer regions of the torus (see the left panel in 
figure H]). The axisymmetric vertical mode, however, shows a more complex behaviour. 
As seen in the right panel of figure HI it involves vertical motions near the inner and 
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Figure 3. The same as figure [TJ but for the axisymmetric radial (left) and vertical 
(right) epicyclic mode in case of a near-extreme Kerr black hole of a — 0.999. 
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Figure 4. Poloidal flow velocity fields of the axisymmetric radial (left) and vertical 
(right) epicyclic mode for a torus of /3 = 0.15 with pressure maximum at ro = 8r g 
orbiting a Kerr black hole with a = 0.5. 



outer edges of the torus, as well as smaller radial flows in regions close to the pressure 
maximum. Its velocity pat tern exhibits the features as on e of the lowest-order slender 
torus modes calculated in iBlaes. Arras and Fragile! (120061 ). the so-called x-mode (see 
their figure 1). The velocity fields of both modes keep their characteristics similar to 
those described above independently of the black hole spin. 



5.2. N on- axisymmetric epicyclic modes 

The lowest-order non-axisymmetric (m = ±1) radial and vertical epicyclic mode 
frequencies are shown in figures [5H3 Again, the (3 = curve in each case refers to 
the test particle frequency and the dash-dotted line specifies the region of allowed 
frequencies (inside the shaded region). Like above, the left panels of figures [5] and [6] 
display the frequencies for tori that orbit a Schwarzschild black hole (a = 0), while 
the corresponding right panels show the same for a rapidly spinning Kerr black hole 
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Figure 5. The non-axisymmctric m — 1 (top ) and m = — 1 (bottom) radial epicyclic 
frequency as a function of r$ for different values of j3 for a M = 1OM black hole. Left 
For a non-rotating black hole. Right For a spinning black hole of a = 0.8. 



(a = 0.8). Figure [7] then displays the frequencies for a near-extreme Kerr black hole 
(a = 0.999). 

The m = 1 radial epicyclic mode frequency of a slightly non-slender torus is given 
by equation ( |52i) . For a test particle (J3 = 0) it reduces to the sum of Qq and the test 
particle axisymmetric radial epicyclic frequency u r Q. For all values of a this frequency 
decreases with increasing torus thickness (see the top panels of figure [5] and the top left 
panel of figure [7]) . 

The m = — 1 radial mode frequency, again given by fl52|) . is in the test particle case 
represented (in absolute value) by the difference between Q and uj r Q. At small radii 
the frequencies for all spin values increase with growing torus thickness (see the bottom 
panels of figure [5] and the bottom left panel of figure [7J. With rising ro, the non-slender 
torus frequencies start oscillating about the (3 = frequency and eventually converge to 
the test particle profile. 

In spherically symmetric spacetimes, the orbital frequency flo and the axisymmetric 
vertical epicyclic frequency ug of a test particle coincide. This is no longer true in case 
of an axially symmetric (Kerr) spacetime or the frequencies of a non-slender torus. The 
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non-axisymmetric m = 1 vertical epicyclic frequency (given by (1561) ) for test particles 
corresponds to the sum of Qq and ujqq. For a non-slender torus the frequency behaves 
up to a < 0.96 similarly to the axisymmetric ojq (compare the top panels of figure [6] to 
figure [2]). Its form slightly differs only for very high spin values (a > 0.96) (see the top 
right panel of figure [7]) . 

The m = — 1 vertical mode frequency (equation fl56|) ) is for a test particle given (in 
absolute value) by the difference between Q and ujqq. In a Schwarzschild spacetime this 
difference equals zero for test particles, but for increased torus thickness the frequency 
grows (bottom left panel of figure [6]). For a Kerr black hole, the frequency converges for 
all spin values to the test particle frequency as the pressure maximum ro moves away 
from the black hole. At smaller radii and for a < 0.96 there is a crossing point with 
the (3 = frequency, such that, in the first interval, rising (3 causes the frequency to 
increase in contrast to the second interval where the frequencies show opposite behaviour 
(figure El bottom right panel). For a > 0.96 the crossing point vanishes and rising torus 
thickness evokes rising frequencies at all radii (bottom right panel of figure Cj). 

Poloidal velocity fields of the non-axisymmetric modes are displayed in figure [8j 
The m = 1 radial mode velocity field exhibits radial flows originating in the central 
regions of the torus and pointing outwards in opposite directions (top left panel). The 
m = — 1 radial mode velocity field has a similar character, except that the radial flows 
point inwards (bottom left panel). The poloidal flow of the m = 1 vertical mode has 
analogous features to that in the axisymmetric case (compare the top right panel of 
figure El to the right panel of figure H]), while for the m — — 1 vertical mode the flow is, 
apart from small variations, mainly vertical (figure [H bottom right panel). Once again, 
the shape of the velocity fields described here remains preserved for all values of the 
black hole spin. 



Blaes et al I (120071 ) presented the properties of the axisymmetric and non- 
axisymmetric radial epicyclic mode frequencies for non-slender tori in a spherical pseudo- 
Newtonian potential, likewise based on calculations accurate to second-order with 
respect to the torus thickness. Comparing their figures 2, 6, 7 to our figures (TJ [2j 
[3j [5] and [7J one may distinguish the behaviour of frequencies calculated in a pseudo- 
Newtonian potential from those calculated in a fully-relativistic Schwarzschild and Kerr 
geometry. Our results show that the radial and vertical mode frequencies in a rotating, 
axially symmetric (Kerr) spacetime follow the same trend as those calculated in a non- 
rotating, spherically symmetric potential, except that the vertical frequency has in a 
rotating spacetime a non-monotonic profile already in the case of test particles. This 
non-monotonicity feature is for thicker tori present in both, axially and spherically 
symmetric, potentials. Generally, only extremely high black hole spin values (a > 0.96) 
cause slight variations in the frequency behaviour. 
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In order to compare the poloidal velocity fields of the axisymmetric modes in Kerr 
ge ometry (our figure [U) to those derived in the pseudo-Newtonian potential see figure 5 



of Blaes et al I (120071 ). 



6. Conclusions 



We have assumed a pressure-supported torus of small radial extent in a Kerr spacetime 
that performs epicyclic oscillations and studied the pressure effects on the epicyclic 
modes properties, i.e., how the modes eigenfunctions and eigenfrequencies of a torus 
differ from those of a free test particle. For this purpose we calculated the relevant 
pressure corrections to the axisymmetric and lowest-order non-axisymmetric epicyclic 
mode eigenfunctions and eigenfrequencies within first- (eigenfunctions) and second-order 
(eigenfrequencies) accuracy in torus thickness. 

In the limit of an infinitely slender torus, the radial and vertical epicyclic oscillations 
occur as global oscillations that correspond to purely radial and vertical displacements 
of the whole torus at epicyclic freq uencies of free test particl e s orbiting at the position 



of th e torus pressure maximum (lAbramowicz et al 



20061 ). Several numerical studies (e.g. iMontero et al 



2006 



2004 



Blaes. Arras and Fragile . 



Rubio-Herrera and Lee 
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Figure 7. The non-axisymmetric m = 1 (top) and m = — 1 (bottom) epicyclic mode 
frequencies in case of a near-extreme Kerr black hole of a = 0.999. Left The radial 
epicyclic mode. Right The vertical epicyclic mode. 



2005 



ISramkova et al 



20071 ) have shown that when the torus gets thicker, its 
(axisymmetric) oscillations occur at lower freque ncies. Th is has been confirmed 



recently by analytic pseudo- Newtonian calculations (IBlaes et al 1 . 120071 ). where pressure 
corrections t o epicyclic modes of a smal l-size, constant specific angular momentum 
torus in the iPaczyhski and Wiital (119801 ) potential were derived. Using the same 
approach, but within the framework of general relativity, we extended their results into 
the rotating Kerr geom etry. For bo t h axi symmetric and the radial non-axisymmetric 
oscillations explored in IBlaes et al I (120071 ). our calculations qualitatively confirm the 
trends as carried out there for a spher ically symmetric potential. As expected (and 
also demonstrated in IBlaes et al I ( 120071 )). the epicyclic mode eigenfunctions of thicker 
tori no longer describe a purely radial or vertical displacement since there appear some 
deviations of the flow in the poloidal velocity fields for both modes. 

The configuration considered here is represented by the idealised model of a 
non-selfgravitating, non-accreting, non-magnetised torus with constant specific angular 
momentumE], studied within a purely hydrodynamical regime. These tori are widely 



+ Uniform specific angular momentum distribution throughout the torus is (aside from simplicity 
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Figure 8. Poloidal velocity fields of the non-axisymmetric epicyclic modes for the 
same torus as in figured] Top The m = 1 radial (left) and vertical (right) epicyclic 
mode. Bottom The m = — 1 radial (left) and vertical (right) epicyclic mode. 



known to be dynami c ally u nstable under global, non-axisymmetric perturbations 
( jPapaloizou fc Pringld . Il984f ) . However, this instab ility can be suppressed when 



1987 



Blaes and Hawleyl . 



accre tion through the inner edge of the torus takes place (IBlaes 
19881 ). A possible interpretation is that torus-like structures can be formed within 
the innermost regions of a standard, nearly-Keplerian accretion disc where several 
physical processes may give rise to pressure gradients that in turn form tori. Such 
torus-like accretion flows seem to appear in three-dimensio nal global MHD simulations 
(foe Villiers and Hawlev et al I . I2OO3I : iMachida etaR hood ). 

Whether purely hydrodynamic modes may exist in the presence of magnetic fields 
that surely play a significant role in the physics of accretion flows is an issue still to 
be inv estigated more deeply . A few studies in that context have been carried out, 
e.g. by iMontero et al I ( 120071 ) who explored the oscillation properties of relativistic tori 
comprising a toroidal magnetic field, and reported similar results as obtained in previous 
investigations of non-magnetised torus oscillations. Their initial set-up, however, did 

reasons) assumed because the relativistic Papaloizou-Pringle equation for n on-constant distributions 
does not describe a self-adjoint eigenvalue problem (see equation (26) in Blaes. Arras and Fragile! 
(2006)), and we do not have the appropriate complete orth onormal set o f eige nfunctions that is necessary 
to apply the perturbation method (see also discussion in Blaes et al I ( 2007 )). 
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not allow for development of the magneto- rotational turbulence (MRI). One of the first 
attempts to investigate the effects of MRI on the pro perties o f hydro dynamic oscillation 
modes via relativistic MHD simulations was done by iFragild (120051 ). 

In general, studying oscillations of black hole accretion discs can improve our 
understanding of the origin of the observed X-ray variability. In particular HF QPOs 
that are detected in the X-ray light curves of several X-ray binaries are often attributed 
to disc oscillations. The results discussed here can be applied to models that directly 
involve epicyclic oscillations. Assuming a particular oscillation model, the identification 
of oscillation frequencies with the frequencies of observed QPO s can provide a p recise 



determination of the mass or spin of the black hol e (e.g. I Wagoner et al I . 12001 



Abramowicz and Kluzniak. 



2001 



Kato and Fukud . 120061 ). Black hole spin estimates for 



several microquasars have been carried out, based on the resonance model considering 
epicyclic osci l lation s in a thin disc that occur at frequencies of free test particles 
(ITorok et al 1 . 120051 ) . Applying pressure corrections to epicyclic fre quencies, our result s 
should be taken into account to obtain (more) accurate estimations (jBlaes et al 1 . 120071 ). 
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Appendix 

Here we display the coefficients introduced to abbreviate the analytic terms we derive 
and use in section HI 



Appendix A. 1. The first- order terms 



The coefficients that appear in the first-order expansion terms of A, Q, f, L in subsection 
12X11 read 
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Appendix A. 2. The second-order terms 



The coefficients in the second-order terms of A, f2, /, L in subsection 14.2.21 take the 
following forms 
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-m 2 (r - 2) - 4mac;f } n } , (A.38) 

L 217 = - L 2W u)ro, (A.39) 
L 218 = -L 21G u 2 d0 . (A.40) 

Appendix A. 3. Eigenf unctions-related coefficients 

The coefficients u^i, w 42 , W51, W52 introduced in tabled] have the form 

^ 2 o(2^o + 2nu; 2 eo - nog 
w 4 i = -2 -2 , (A.41) 

g|o(2g|o + 2n^ r 2 - nof) 
W42 - -2 — 2 , (A.42J 

^(90 ~~ ^rO 

Ufc = _^0( 2 ^0_+ 2 ^f0-^5 2 ) ; (A43) 
w 52 - -2 -2 • (A.44) 
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